Non-linear macroscopic polarization in III-V nitride alloys 
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We study the dependence of macroscopic polarization on composition and strain in wurtzite III-V 
nitride ternary alloys using ab initio density-functional techniques. The spontaneous polarization 
is characterized by a large bowing, strongly dependent on the alloy microscopic structure. The 
bowing is due to the different response of the bulk binaries to hydrostatic pressure, and to internal 
strain effects (bond alternation). Disorder effects are instead minor. Deviations from parabolicity 
(simple bowing) are of order 10 % in the most extreme case of AlInN alloy, much less at all other 
compositions. Piezoelectric polarization is also strongly non-linear. At variance with the spontaneous 
component, this behavior is independent of microscopic alloy structure or disorder effects, and due 
entirely to the non-linear strain dependence of the bulk piezoelectric response. It is thus possible 
to predict the piezoelectric polarization for any alloy composition using the piezoelectricity of the 
parent binaries. 
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I. INTRODUCTION 

III-V nitrides and their alloys have recently emerged 
as a strategic materials system for the design and manu- 
facture of high-freauency light emitting devices and high, 
power transistors.El Recent theoretical investigationsBcl 
have suggested that the macroscopic polarization of 
wurtzite nitrides has a major influence on the design 
criteria of low-dimensional nanostructures aimed to ad- 
dress the technological goals mentioned above. Indeed, 
III-V nitrides in their wurtzite structure posses a large 
spontaneous polarization and larger-than-usual piezo- 
electric constants.0 As a consequence, fixed and approx- 
imately two-dimensional charges resulting from polar^. 
ization discontinuities occur at nitride heterointerfacesjj 
causing built-in electrostatic fields and/or free-charge 
accumulation, completely unusual in conventional III- 
V's.tJo By andJarge, these findings have been confirmed 
experimentally!] through the observation of optical shifts 
in quantum wells, and of high-density two-dimensional 
electron gases in field-effect transistor structures, and 
have been used as a guideline in the design and optimiza- 
tion of non-conventional high-frequency, high-breakdown 
voltage heterostructure field-effect transistors. 

In nanostructure modeling or experiment interpreta- 
tion, it was assumed so farta that polarization in nitride 
alloys would interpolate linearly between the limiting, val- 
ues determined by the parent binary compounds J3 i.e. 
follow a Vegard-like behavior. The structural properties 
generally follow Vegard's law in many semiconductor al- 
loys. On the other hand, electronic properties such as the 
optical band gap often depart from linear behavior: it is 
thus a plausible expectation that macroscopic polariza- 
tion may also exhibit non-linearity. In this paper, we in- 
vestigate the dependence of the macroscopic polarization 
on the composition and microscopic structure in AlGaN, 
InGaN and AlInN alloys, by means of ab initio density- 
functional techniques. Whereas the structural properties 



do follow Vegard's law closely, wc find a strong non-linear 
dependence of the polarization on alloy composition. The 
non-linearity of the spontaneous polarization depends ap- 
preciably on the microscopic structure of the alloy. By 
contrast, the non-linearity of the piezoelectric polariza- 
tion is a pure bulk effect: Vegard's law holds closely if 
the non-linear strain dependence of bulk piezoelectricity 
is accounted for. 



II. ALLOY STRUCTURE 

The calculation of the macroscopic polarization of an 
alloy requires the knowledge of its equilibrium structure. 
Nitride alloys are commonly thought to be characterized 
by a random distribution of the group-Ill elements on the 
wurtzite cation sites, while anion sites-.are always occu- 
pied by nitrogen. Regent experimental^ work (with sup- 
port from theoreticafl investigations) has suggested that 
ordering takes place in some nitride alloys under suit- 
able growth conditions. The ordering in question was re- 
ported to consist of 1 x 1 super-lattice structures oriented 
along the (0001) direction, -aath Al (or In) atomic planes 
alternating with Ga planes.Q Given these suggestions, a 
theoretical investigation of the alloy phase diagram would 
require the study of surface thermodynamics and kinet- 
ics as a function of the growth conditions. This is out- 
side the scope of the present work. Here we accept as 
plausible that both random and ordered structures may 
be obtained depending on the growth conditions, and 
consider a sample of both structures in our polarization 
calculation. The ordered structure just mentioned bears 
evident similarities with the CuPt-like structure found in 
epitaxially ordered alloys of zincblende materials,!] where 
the super-lattices are formed in the (111) direction; we 
therefore label this structure as CP. On the other hand, 
the microscopic structure of a random alloy can also be 
represented efficiently in periodic conditions using the 
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special quasi-random structure (SQS) method.tl Given its 
successes in the prediction of band gap non-linearity in 
zincblende-structure alloys, the SQS method is the nat- 
ural choice to study polarization non-linearity in III-V 
nitride alloys. We enforce the periodic boundary condi- 
tions needed to compute thp macroscopic polarization in 
the Berry phase approachL2l by using repeated supercells, 
as done previously in Ref. |ll| to study Born charges and 
IR spectra in amorphous SiC>2. As a compromise be- 
tween computational workload and the accurate descrip- 
tion of random structures, a 32-atom 2x2x2 wurtzite 
supercell is adopted. In analogy to zincblende-based al- 
loys, it is possible to reproduce the statistical properties 
(pair correlation functions, etc.) of a random wurtzite 
alloy for a molar fraction x=0.5 by suitably placing the 
cations on the 16 sites available in the cell, i.e. by a so- 
called SQS-16 structure. Eves, smaller structures such as 
SQS-8 were previously founcLJ to be adequate to mimic 
random zincblende-based alloys, so we are confident that 
our SQS-16 is a good model of the present random al- 
loys. Since other molar fractions cannot be described as 
easily, and also because non-linear effects are expected 
to be largest for this concentration, our investigation on 
random alloys is restricted to the case £=0.5. 

As a -further ingredient, we consider ordered 
structureaj other than the CuPt mentioned earlier for 
each of the ternaries (AlGaN, InGaN, AlInN). The 
chalchopyrite-like structure (henceforth labeled CH) is 
defined by each anion site being surrounded by two 
cations of one species and two of the other, with the con- 
straint of fitting periodically into our (2x2x2) wurtzite 
supercell. This structure is highly symmetric, as there 
are only two kinds of inequivalent anion sites, differing 
in the orientation of the neighbors, and not in their chem- 
ical identity. Among the possible ordered structure, CH 
is in a sense the most homogeneous for the given compo- 
sition. A further useful ordered structure we consider is a 
luzonite-like structure, resembling the luzonite structured 
used for zincblende-based alloys. In this structure (la- 
beled LZ), each nitrogen atoms is surrounded by three 
cations of one species, and one of the other: in a sense, 
this is the analog to the CH structure for molar fractions 
x = 0.25 and 0.75. 

The comparison between chalchopyrite-like and SQS 
"random" structures will give us insights on the effect 
of randomness versus ordering without the biases due to 
specific super-lattice ordering as in the CuPt structure. 
Luzonite-like structure will provide values of the polar- 
ization at intermediate molar fractions. 



III. METHOD 

As for the technical ingredients of the calculations, 
we work within the density-functional-theory pseudopo- 
tential planerwave framework as implemented in the 
VASP codc.EJ Ultrasoft pseudopotentialala provided with 



vASpHa are used for all the elements involved; Ga and In 
d semicore states are included in the valence. Since we 
are interested in the effects on polarization due to (among 
other) internal strains related to size mismatch of the 
two alloyed materials, it is imperative that our calcula- 
tion reproduce relative mismatches with the highest pos- 
sible accuracy. Comparing DFT calculations using both 
the LDA and GGA (in its PW91 variantPJ) exchange- 
correlation functional, we found that GGA does a bet- 
ter job of reproducing the relative mismatch between the 
binary constituents, besides getting closer to the exper- 
imental structure. ta We therefore used the GGA for all 
the present calculations. We reported recently in some 
detail the properties of binary nitrides, including polar- 
ization, as calculated in the GGA vs LDA in Refs. and 

Another precaution adopted throughout this work was 
to use the same 32-sites supercell for both binaries and al- 
loys, the same k-points for reciprocal space integration (a 
4x4x4 Monkhorst-Pack mesh) and for the Berry phase 
calculations, and the same safe plane-wave cutoff of 325 
eV, regardless of composition. Although more costly, 
this route reduces the systematic convergence errors. All 
forces and stress components have been converged to zero 
within 0.005 eV/A and 0.01 Kbar, respectively, for free- 
standing alloys. For GaN-epitaxial alloys (Sec. 0), only 
the zz stress component has been zeroed. The polariza- 
tion of the binary compounds as obtained in the supercell 
calculations is always within less than 0.001 C/m 2 of the 
value obtained in a bulk calculation using a 12x12x12 
k-points mesh. 

IV. SPONTANEOUS POLARIZATION 

We start our analysis comparing the values of the lat- 
tice constants and spontaneous polarization computed 
for the free standing alloys, with the predictions of Ve- 
gard's law. The latter relates the lattice constants a 
and c, and the spontaneous polarization P sp of an al- 
loy A x Bi_ x N, to their values for the binary costituents 
AN and BN (A and B indicate two generic cations among 
In, Al, and Ga) by 

a (A x Bi_ x N) = x a (AN) + (1 - x) a(BN), (1) 

c (A x Bi_ x N) = x c(AN) + (1 - x) c(BN), (2) 

and 

P S p(A x B!_ x N) = xP sp (AN) + (1 -x)P sp (BN). (3) 

The existence of compositional non-linearities (i.e. non- 
Vegard behavior) in a generic measurable quantity in 
an alloy can be described in a first approximation by 
a parabolic model involving a so-called bowing parame- 
ter. In particular, we expect to be able to describe the 
possible non- linear dependence of the spontaneous polar- 
ization on the composition as 
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P(A x Bi_ x N) = xP(AN) + (1 - x) P(BN) 

-6 A Bx(l-a;), (4) 

where the bowing parameter is by definition as 

6ab = 2P(AN) + 2P(BN) - 4P(A .5oBo. 50 N) (5) 

involving only the knowledge of the polarization for 
ternary alloys with molar fraction x — 0.5. Given the 
bowing parameter, the spontaneous polarization can be 
obtained at any composition. 



A. Polarization non-linearity 

In Fig. [l], we compare with Vegard-law predictions our 
calculated equilibrium lattice parameters a and c for all 
the alloy structures and compositions studied. Both a 
and c follow quite closely Vegard's law in all cases. The 
dependence of the polarization on composition is then 
the same as that on the lattice parametcr(s), modulo a 
multiplicative factor. Therefore, in Figs.|| and || we can 
compare with Vegard-law predictions the calculated data, 
respectively, for the random and Cu-Pt-ordered (FigJ|), 
and chalcopyrite- or luzonite-like ordered (Fig. ||) alloys, 
displayed in the (P,a) plane. The ordinates of the filled 
symbols are the calculated spontaneous polarization, the 
abscissas are the calculated lattice constants; the open 
symbols and dashed lines are the Vegard-law predictions; 
the solid lines are fits of Eq.[| to the calculated values. 
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FIG. 1. Basal (a) and axial (c) lattice constants of 
wurtzite-based nitride alloys. Directly calculated values for 
various alloy structures are denoted by open circles (random) , 
squares (CH and LZ), triangles (CP). The dashed lines are 
Vegard's law. 



We identify two main points. First, there is a large 
non-linear dependence of the spontaneous polarization 
on the lattice parameter, hence on alloy composition. Al- 
loys with large lattice mismatch between the constituents 
(e.g. AlInN and InGaN) show the largest bowing, while 
in AIGaN the non-linearity is modest, although sizable. 
Regardless of the alloy structure and composition, the 
spontaneous polarization shows always an upward bow- 
ing. Second, CP-ordered (Fig. ||) and CH-ordered (Fig. 
^) alloys have very different bowing. This is at variance 
to the case of the band gap: ordering is known to cause 
an increase of the band gap bowing, but more homoge- 
neous across different ordered structures, and never as 
dramatic as that of polarization in the CP ordered alloy. 
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FIG. 2. Spontaneous polarization versus equilibrium lat- 
tice constant in free-standing random and CuPt-ordered ni- 
tride alloys (filled symbols). Solid lines are interpolations ac- 
cording to Eq. ^. Dashed lines and open symbols are Vegard 
predictions. 

The CH and LZ ordered structures offer the possibility 
to test the validity of the interpolation model based on 
the bowing parameter. The ordered LZ structure is anal- 
ogous to CH for molar fractions of 0.25 and 0.75. The size 
of the deviations of the LZ values from those predicted by 
the fit of Eq|I| through the binaries and the CH (x=0.5) 
value, indicates whether or not non-parabolicity occurs 
in the P(x) relation for CH-like order. Since the polar- 
ization of the CH structure behaves qualitatively as that 
of the random structure (Fig. |[), the conclusions for CH 
may be exported to the random phase. In Fig. [|the val- 
ues of the polarization calculated for the LZ structures 
at molar fractions 0.25 and 0.75 lay very close to the 
parabolic curve for InGaN and AIGaN; this confirms the 
accuracy of the quadratic hypothesis embodied in Eq^. 
For AlInN, the calculated values lay somewhat above (be- 
low) the parabolic relation for £=0.25 (x=0.75), show- 
ing that some polarization non-parabolicity occurs, and 
specifically that the bowing is higher for low In concentra- 
tion in AlInN. Compared to the quadratic non-linearity 
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(Eq.[|) , the non-parabolicity causes relatively modest ad- 
ditional deviations (of order 10 %) for AlInN. A legiti- 
mate conclusion is then that the bowing relation should 
predict polarization in AlInN with about 10 % accuracy, 
and much better for InGaN and AlGaN. 
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FIG. 3. Spontaneous polarization versus equilibrium lat- 
tice constant for free-standing chalcopyrite- and luzonite-like 
ordered nitride alloys (filled symbols; open symbols are used 
for InGaN for clarity) . Solid lines are interpolations according 
to Eq. Dashed lines are Vegard predictions. The shaded 
triangle represents the polarization for the random alloys as 
reported in Fig. ^. 

It is clear that these result have a considerable tech- 
nological relevance for nanostructure devices. Polariza- 
tion differences across interfaces in, say, naulti-quantum- 
wells, generate localized interface chargesH which in turn 
cause electrostatic fields and/or charge accumulation at 
the interfaces.Qti The non-linearity just reported modifies 
considerably the polarization-induced interface charge 
densities, and must be accounted for in any attempt at 
quantitative comparisons with, or predictions on exper- 
iments. We will be considering this as pect in detail in 
forthcoming work. In the following Sec. IV B we attempt 



instead to gain insight on the nature of the non-linear 
behavior of the polarization in III-V nitrides. 



B. Microscopic origin of polarization non-linearity 

It is known from eaelier theoretical work on III-V 
zincblende-based alloyst-3 that band gap non-linearities 
can be attributed to essentially three sources: i) chemi- 
cal effects due to the different cation electronegativity of 
the components; ii) internal strain effects, due to varying 
cation-anion bond lengths (the so-called bond alterna- 
tion); Hi) disorder effects due to the random distribution 
of the chemical elements on the cation sites. Chemical 
effects and internal strain are dominant, whereas disor- 
der gives only a modest contribution. On the basis of 
these indications, one expects a large bowing in random 



alloys whose constituents have a large lattice mismatch. 
Such a behavior for the fundamental gap bowing,was in- 
deed found for zincblende III-V nitrides alloysta Also, 
the bowing in ordered structures should be similar to, 
but larger than in random structures. 

For the polarization, we only find a qualitative behav- 
ior similar to the band gap bowing in the case of random 
alloys, while the ordered phases CH and CP do not seem 
to show identifiable trends. It is especially surprising that 
i) the random and CP phases of AlGaN have nearly the 
same bowing; ii) the InGaN CH ordered alloy shows a 
smaller bowing than its random phase. The reason for 
this anomalous behavior must be found in the peculiar 
nature of the polarization bowing. While the band gap 
is a scalar, the polarization is a vector of fixed direction, 
namely, for wurtzites, the (0001) direction. Thus, bond 
alternation will affect the polarization bowing only if it 
changes the projection of the bond length in the (0001) 
direction. This idea is supported by the fact that in pure 
binaries the polarization is strongly influenced by the rel- 
ative displacement of-the cation and anion sub-lattices in 
the (0001) direction.!^ 

Indeed, we find that a clear correlation exists between 
the so-called u parameter of the wurtzite structure, i.e. 
the bond length along the singular polar (or pyroelectric) 
axis, and the value of the polarization. In Fig. ^, we show 
the calculated spontaneous polarization of free-standing 
binaries and x=0.5-alloys in the various structures, vs. 
the average internal parameter u. The latter is defined 
as the average value of the projection of the connecting 
vector of a nitrogen atom with its first neighbor in the 
(0001) direction along this same direction!^ This defini- 
tion can be used also for random phase alloys in spite of 
the displacements off the ideal sites. We see in Fig. ^ 
that for a given alloy composition, the spontaneous po- 
larization of free-standing alloys of different microscopic 
structure depends linearly on the average u parameter 
of the alloy structure. This indicates that spontaneous 
polarization differences between alloys of the same com- 
position are due essentially to structural and bond al- 
ternation effects; disorder appears to have a negligible 
influence. It appear instead that the chemical identity 
of the constituents plays a role of some importance (be- 
yond the obvious fact that it implicitly determines the 
structure). 

The points mentioned above are now easily explained. 
The random and CP phases in AlGaN have almost the 
same average u, hence almost the same polarization. In 
InGaN, the random alloy has a larger u than the CH 
phase, while the opposite holds for AlInN and AlGaN: 
therefore the CH-vs-random bowing behavior in InGaN is 
opposite with respect to the others. Also, the huge bow- 
ing of CP-ordered AlInN and InGaN matches the very 
large deviation of the average u in those phases from the 
random and CH-like structures. If structure were the 
only source of polarization bowing, all of the points in 
Fig. H would fall on the same straight line. 
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FIG. 4. Spontaneous polarization versus average u param- 
eter in ternary alloys. Open circles, squares, triangles, refer 
to random, CH-like, CP-like structures, respectively. 

To check this role, we set up a model based on the 
polarization of binaries and alloys in a constrained ideal 
wurtzite structure. In this structure, only the a parame- 
ter is independent, whereas c and u are fixed at the values 
determined by maximal sphere packing, namely w=0.375 
and c/a—^J%J?>. Each nitrogen atom is then surrounded 
by four equidistant cations - that is, all bonds have the 
same length for a given lattice constant a. By construc- 
tion, then, bond alternation do not play any role, and 
the effects of chemical identity of the constituents can be 
easily disentangled. 

We assume that Vegard's law holds for the lattice con- 
stant a, as was found to be the case for the alloys in 
the preceeding Subsection. This establishes a (linear) 
relation between composition and lattice constant. We 
then calculate the polarization in each of the binary ni- 
trides in their ideal structure as function of the lattice 
constant a(x). Finally, we express the alloy polarization 
as a composition- weighted Vegard-like average of the po- 
larizations of the binaries as 

P sp (A x B^ x N) = xP^(AN) + (1 - x)P< x \BN). (6) 



In this model, any non-linearity must come from the dif- 
ferent response of polarization to changes in a(x), hence 
to hydrostatic compression, in the ideal binary com- 
pounds. 

We report in Fig. [s] three sets of quantities: the po- 
larizations computed in the ideal wurtzite structure for 
the three ternary alloys at their calculated equilibrium 
volume (filled circles); the polarization of the binaries in 
the ideal structure as a function of the lattice parameter 
a (open symbols); the Vegard interpolation EqJ^ of the 
latter polarizations (the bowed shaded triangle) for the 
alloys. We see that the computed values for the alloys and 
the Vegard prediction essentially coincide. Therefore, the 
origin of the "chemical" non-linearity, and its large val- 
ues in In-containing alloys, becomes clear: in A1N and 
GaN, the polarization decreases with hydrostatic pres- 



sure, while it increases in InN. Also to be noted, polar- 
izations in the ideal structure are between 35 to 50 % of 
their value in free standing alloys, and despite the ab- 
sence of bond alternation, the bowing is still very large. 
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Solid lines are Vegard interpolations of the latter values. This 
shows that the effect of hydrostatic pressure on bulk polar- 
izations reflects directly on the polarization in the alloys. 

The above model provides an expression of the bowing 
parameter of ideal wurtzite structure alloys as function 
of the polarization response to hydrostatic pressure: 
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The agreement of the latter expression with the b result- 
ing from a fit to the calculated values is very good (e.g. 
for the extreme case of AlInN, 6 mo doi=-0.0225 C/m 2 , 
while from direct calculation we get -0.0208 C/m/ 2 ). On 
the basis of the model, it is now understandable that the 
AIGaN bowing is pretty moderate since the region of in- 
terest is small (3.1 A to 3.2 A) and the responses to 
hydrostatic pressure of AIN and GaN are similar. On 
the contrary, in the large range 3.1-3.6 A, AIN and InN 
have opposite behaviors, whence the huge bowing found 
in AlInN alloys. The same goes, although to a lesser 
extent, for InGaN alloys. 



V. PIEZOELECTRIC POLARIZATION 



The knowledge of the spontaneous polarization in al- 
loys is not sufficient to predict the values of interface 
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charge accumulations or electrostatic fields in wurtzite- 
nitride-based nanostructures, because the latter are usu- 
ally grown pseudomorphically and under strain on the 
(0001) or (0001) faces of a substrate (typically sapphire 
or SiC, with or without A1N or GaN buffer layers). The 
ensuing symmetry-conserving strain causes a change in 
polarization that amounts to a piezoelectric polarization. 
In general, the polarization in the active layers of a nanos- 
tructure is a combination of spontaneous and piezoelec- 
tric polarization. In this context, the aim of this Section 
is threefold. First, we demonstrate that piezoelectricity 
in nitride alloys is non-linear; second, we show that this 
non-linearity is due to a pure bulk effect, i.e. the non- 
linear behavior of bulk binary piezoelectric constants vs 
symmetry-conserving strain; third, we suggest how to use 
this understanding in practice to predict piezoelectricity 
in any nitride alloy. 

We consider for definitcness (but with no loss of gen- 
erality as to the scheme to predict piezoelectricity in 
arbitrary alloys) the technologically relevant case of an 
alloy pseudomorphically grown on an unstrained GaN 
substrate. Using the same technical ingredients as pre- 
viously, we repeat the polarization calculations with the 
costraint that a alloy =a GaN , fully reoptimizing all struc- 
tures. We compute the piezo component as the difference 
of the total polarization obtained in this calculation and 
the spontaneous part calculated previously for the free- 
standing alloy. 

In Fig. [| we show the piezoelectric polarization as a 
function of the alloy composition. Symbols are the cal- 
culated polarizations in the various alloy structures and 
compositions; it is clear that, contrary to spontaneous 
part, the piezoelectric polarization component hardly de- 
pends on the microscopic structure of the alloy. We now 
ask whether the piezoelectric polarization of the alloy can 
be reproduced by a Vegard-like model containing only the 
knowledge of the properties of the binaries, i.e. 
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FIG. 6. Piezoelectric component of the macroscopic po- 
larization in ternary alloys epitaxially strained on a GaN 
substrate. Open symbols are directly calculated values for 
random (circles), CH-like and LZ-like (squares), and CP-like 
(triangles) structures, respectively. Dashed lines are the pre- 
diction of linear piezoelectricity, while the solid lines is the 
prediction of EqiS using the non-linear bulk polarization as 
reported in Fig. m. 

This is due to a genuine non-linearity of the bulk piezo- 
electricity of the binary constituents, of non-structural 
origin; hence, the bowing due to the microscopic struc- 
ture of the alloys is negligible. The proof of this state- 
ment is in two steps. First, we directly calculate the 
piezoelectric polarization as a function of basal strain for 
AIN, GaN, and InN, optimizing all structural parame- 
ters. The results, collected in Fig. ^, clearly show that 
the piezoelectric polarization of the binaries is an ap- 
preciably non-linear function of a, i.e. of basal strain. 
Since all lattice parameters closely follow Vegard's law, 
the non-linearity cannot be related to deviations from 
linearity in the structure. 
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BN 
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In a first approximation, one may calculate the piezoelec- 
tric polarization of the binary compounds for symmetry- 
conserving in-plane and axial strains as 



2e 33 ei, 



(9) 



where the piezoelectric constants e are calculatedo in the 
equilibrium structure of the binary AN, and by definition 
do not depend on strain. The dashed lines in Fig. [| rep- 
resent the piezoelectric term as computed from the above 
relations using the piezoelectric constants computedEI for 
the binaries. Clearly, when combined with Eq!fy Vegard's 
law (Eqj8j) fails to reproduce the calculated polarization, 
and misses the strong non-linearity of the piezoelectric 
term evident in Fig. ffl. 
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FIG. 7. Piezoelectric polarization in binary nitrides as 
function of basal strain (symbols and solid lines) compared 
to linear piezoelectricity prediction (dashed lines) . The c and 
u lattice parameters are optimized for each strain. 



G 



Second, we plug the non-linear piezoelectric polariza- 
tion just computed for the binaries into the Vegard inter- 
polation Eq|| We obtain thereby the solid lines in Fig. 
|6f we see that the predictions of this scheme are in excel- 
lent agreement with the polarization calculated directly 
for the alloys. Therefore, we can conclude that a) the 
non-linearity in bulk piezoelectricity dominates over any 
effects related to disorder, structure, bond alternation, 
etc., and that b) Vegard's law holds for the piezoelec- 
tric polarization of III-V nitrides alloys, provided that the 
non-linearity of the bulk piezoelectric of the constituents 
is accounted for. 

The latter conclusion suggests a straightforward 
scheme for predicting the piezoelectric polarization of any 
nitride alloy at any strain: pick a value for x, calculate 
the basal strain e(x) from Vegard's law, calculate P pz 
from EqJ| using the non-linear piezoelectric polarization 
of the binaries (Fig. |?]). This scheme is of obvious inter- 
est in the modeling of nanostructures, especially for high 
In contents and AlInN alloys. 

VI. SUMMARY AND ACKNOWLEDGEMENTS 

In this paper, we studied the non-linear dependence 
of the polarization on composition and strain in ternary 
III-V nitride wurtzite-based alloys. The spontaneous po- 
larization is characterized by a large bowing, strongly 
dependent on the alloy microscopic structure. The bow- 
ing is due to the different response of the bulk binaries to 
hydrostatic pressure, and to internal strain effects due to 
bond alternation. Disorder effects are minor. Deviations 
from parabolicity (simple bowing) are of order 10 % in 
the most extreme case (the AlInN alloy), much less for 
all other compositions. 

Piezoelectric polarization is also strongly non-linear. 
This behavior is independent of microscopic structure or 
disorder effects (the structural parameters for all alloys 
closely follow Vegard's law and linear elasticity), and is 
due entirely to the non-linear strain dependence of the 
bulk piezoelectric response. Based on our understand- 
ing of this behavior, we suggested that it is possible to 
predict the (non-linear) piezoelectric polarization for any 
alloy composition using the knowledge of the binary com- 
pound piezoelectricity. 

We acknowledge support from a MURST Cofin99 
project and from Iniziativa Trasversale Calcolo Parallelo 
of INFM. 
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